Quantum-Classical Dynamics of Wave Fields 
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An approach to the quantum-classical mechanics of phase space dependent operators, which has 
been proposed recently, is remodeled as a formalism for wave fields. Such wave fields obey a system 
of coupled non-linear equations that can be written by means of a suitable non-Hamiltonian bracket. 
As an example, the theory is applied to the relaxation dynamics of the spin-boson model. In the 
adiabatic limit, a good agreement with calculations performed by the operator approach is obtained. 
Moreover, the theory proposed in this paper can take nonadiabatic effects into account without 
resorting to surface-hopping approximations. Hence, the results obtained follow qualitatively those 
of previous surface-hopping calculations and increase by a factor of (at least) two the time length 
over which nonadiabatic dynamics can be propagated with small statistical errors. Moreover, it is 
worth to note that the dynamics of quantum-classical wave fields here proposed is a straightforward 
non-Hamiltonian generalization of the formalism for non-linear quantum mechanics that Weinberg 
introduced recently. 



I. INTRODUCTION 

There are many instances where a quantum-classical 
description can be a useful approximation to full quan- 
tum dynamics. Typically, a quantum-classical picture 
often allows one to implement calculable algorithms on 
computers whenever charge transfer is considered within 
complex environments, such as those provided by pro- 
teins or nano-systems in general With respect 
to this, an algebraic approach has been recently pro- 
posed [1, Q in order to formulate the dynamics and 
the statistical mechanics [3| of quantum-classical sys- 
tems. General questions regarding the quantum-classical 
correspondence have also been addressed within a sim- 
ilar framework [5]. The approach of Refs. [3, Q rep- 
resents quantum-classical dynamics by means of suit- 
able brackets of phase space dependent operators and 
describes consistently the back-reaction between quan- 
tum and classical degrees of freedom. Notably, a partic- 
ular implementation of this formalism has been used to 
calculate nonadiabatic rate constants in systems model- 
ing chemical reactions in the condensed phase How- 
ever, such schemes have only permitted the simulation of 
short-time nonadiabatic dynamics because of the time- 
growing statistical error of the algorithm. Neverthe- 
less, the algebraic approach 0, [1], underlying the algo- 
rithms of Refs. [gJ, has some very nice features, such as 
the (above mentioned) proper description of the back- 
reaction between degrees of freedom, that one should 
not give up when addressing quantum-classical statis- 
tical mechanics. Moreover, quantum-classical brackets 
define a non-Hamiltonian algebra |3] so that their ma- 
trix structure allows one to introduce quantum-classical 
Nose-Hoover dynamics Q and to define the statistical 
mechanics of quantum-classical systems with holonomic 
constraints ^ . All of the above features of the formalism 
are highly desirable when studying complex systems in 



*E-mail: ascrgi@unime.it 



condensed phases. Therefore, it is worth to search for a 
reformulation of the theory of Refs. [1, 0, 11] that, while 
maintaining such features, could be used to integrate re- 
liably long-time nonadiabatic dynamics. 

To this end, one can note that, within standard quan- 
tum mechanics, some problems that are formidable to 
solve by means of the dynamics of operators become 
much simpler to handle when, instead, the time evolu- 
tion of wave functions is considered Q . Hence, for anal- 
ogy, it might also happen that, within quantum-classical 
mechanics, the correspondence between operators and 
quantum-classical wave functions could open new possi- 
bilities for useful approximations in order to carry long- 
time calculations efficiently. Indeed, finding and applying 
the correspondence between operator and wave scheme of 
motion in quantum-classical mechanics is the scope of the 
present paper. A wave picture for quantum-classical dy- 
namics can be found by direct algebraic manipulation of 
the equation of motion for the density matrix. In prac- 
tice, the single equation obeyed by the quantum-classical 
density matrix is mapped onto two coupled non-linear 
equations for quantum-classical wave fields. Despite its 
non-linear character, such a quantum-classical dynamics 
of phase space dependent wave fields corresponds exactly 
to the dynamics of phase space dependent operators dis- 
cussed in Refs. 0, H, 0, 0, ll| and can be used to devise 
novel algorithms and approximation schemes. 

The abstract algebraic equations here presented are 
readily expressed in the adiabatic basis and applied, in 
order to provide an illustrative example, to the spin- 
boson model and its relaxation dynamics both in the 
adiabatic and nonadiabatic limit. By making a suit- 
able equilibrium approximation to the non-linear wave 
equations, it is found that nonadiabatic dynamics can be 
propagated, within the wave picture, for time intervals 
that are a factor of two-three longer than those which 
have been spanned in Ref. by means of the operator 
theory @,i,i0,i. Such a result is very encouraging 
for pursuing the long-time integration of the nonadiabatic 
dynamics of complex systems in condensed phases. 

Following a line of research that investigates the rela- 
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tions between classical and quantum theories [llj, it is 
worth to note that the wave picture of quantum-classical 
mechanics, which is introduced in this paper, generalizes 
within a non-Hamiltonian framework the elegant formal- 
ism that Weinberg [12] proposed for describing possible 
non-linear effects in quantum mechanics 13]- 

This paper is organized as follows. In Section HIl the 
non-Hamiltonian algebra of phase space dependent oper- 
ators is briefly summarized. In Section Hill the quantum- 
classical dynamics of operators is transformed into a the- 
ory for phase space dependent wave fields evolving in 
time. Such a theory for wave fields is also expressed 
by means of suitable non-Hamiltonian brackets: in this 
way a link is found with the generalization of Weinberg's 
non-linear formalism given in Appendix [Xj More specif- 
ically, in Appendix [XJ Weinberg's formalism is briefly 
reviewed and its symplectic structure is unveiled. Then, 
this structure is generalized by means of non-Hamiltonian 
brackets. Therefore, one can appreciate how the gen- 
eralized Weinberg's formalism establishes a more com- 
prehensive mathematical framework for non-linear equa- 
tions of motion, comprising phase space dependent wave 
fields as a special case. In Section IIVI the abstract non- 
linear equations of motion for quantum-classical fields 
are represented in the adiabatic basis and some consider- 
ations, which pertain to the numerical implementation, 
are made. By making an equilibrium ansatz, in SectionlVl 
the non-linear equations of motion are put into a linear 
form and the theory is applied to the spin-boson model. 
Section is devoted to conclusions and perspectives. 



II. NON-HAMILTONIAN MECHANICS OF 
QUANTUM-CLASSICAL OPERATORS 

A quantum-classical system is composed of both quan- 
tum X and classical X degrees of freedom, where X = 
(i?, P) is the phase space point, with R and P coordi- 
nates and momenta, respectively. Within the operator 
formalism of Refs. [1, H 0, H], the quantum variables 
depend from the classical point, X, of phase space. The 
energy of the system is defined in terms of a Hamiltonian 
operator H = H{X), which couples quantum and clas- 
sical variables, by £' = Ti' J dX H{X). The dynamical 
evolution of a quantum-classical operator x{X) is given 
by Si 



where 



x{X,t) = ^[H,x{X,t)\^-^[H,x{X,t)} 



+ ^[x{X,t),H}^=(^H,x{X,t)) , (1) 



H,X 



B 



H 

X 



(2) 



is the commutator and 
{H,x}b = 



^dH^ dx_ 

2^ dxr^dx, 



(3) 



is the Poisson bracket [l^ . Both the commutator and the 
Poisson bracket are defined in terms of the antisymmetric 
matrix 



1 
-1 



(4) 



The last equality in Eq. ([T]) defines the quantum-classical 
bracket. Following Refs. 7, 8, 15], the quantum-classical 
law of motion can be easily casted in matrix form as 



d i r - n 



dt 
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where 
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(5) 



(6) 



The structure of Eq. ([5]) is that of a non-Hamiltonian 
commutator, which will be defined below in Eq. (jlOp . 
and as such generalizes the standard quantum law of mo- 
tion [7,]. The antisymmetric super-operator in Eq. ([5]) 
introduces a novel mathematical structure that charac- 
terizes the time evolution of quantum-classical systems. 
The Jacobi relation in quantum-classical dynamics is 
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X, 
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The explicit expression of has been given in Ref. 0] 
where it was shown that it may be different from zero 
at least in some point X of phase space: for this reason 
the quantum-classical theory of Refs. [1, H, 0, 9 can be 
classified as a non-Hamiltonian theory. 

It is worth to note that the quantum-classical law of 
motion in Eq. ([5]) is a particular example of a more gen- 
eral form of quantum mechanics where time evolution is 
defined by means of non-Hamiltonian commutators. The 
non-Hamiltonian commutator between two arbitrary op- 
erators X a-iid ^ is defined by 



[x,iu^[x i]-^ 



(8) 



where is an antisymmetric matrix operator of the form 

-m] 



(9) 



where /[?)] can be another arbitrary operator or func- 
tional of operators. Then, generalized equations of mo- 
tion can be defined as 



= i{H,x]n. 



H 

X 



(10) 
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The non-Hamiltonian commutator of Eq. ([8]) defines a 
generalized form of quantum mechanics where, neverthe- 
less, the Hamiltonian operator H is still a constant of 
motion because of the antisymmetry of Jl. 



III. QUANTUM-CLASSICAL WAVE 
DYNAMICS 

In Refs. 0, 0], quantum-classical evolution has been 
formulated in terms of phase space dependent operators. 
In this scheme of motion operators evolve according to 



x{X,t) = expU 



exp{ii£}x(X) , 



(11) 



where the last equality defines the quantum-classical Li- 
ouville propagator. Quantum-classical averages are cal- 
culated as 



(X)(0 = Tr'y dXp{X)x{X,t) 

= Tr' / dXp{X,t)x{X) , (12) 



where p{X^ is the quantum-classical density matrix and 
p{X,t) = exp{—itC}p{X). Either evolving the dynam- 
ical variables or the density matrix, one is still dealing 
with phase space dependent operators: viz., one deals 
with a form of generalized quantum-classical matrix me- 
chanics. As it has been discussed in the Introduction, 
this theory has interesting formal features and a cer- 
tain number of numerical schemes have been proposed 
to integrate the dynamics and calculate correlation func- 
tions However, these algorithms have been ap- 
plied with success only to short-time dynamics because of 
statistical uncertainties that grow with time beyond nu- 
merical tolerance. With this in mind, it is interesting to 
see which features are found when the quantum-classical 
theory of Refs. 0, Q is mapped onto a scheme of mo- 
tion where phase space dependent wave fields, instead of 
operators, are used to represent the dynamics. 

As it is well known in standard quantum mechan- 
ics, the correspondence between dynamics in the Heisen- 
berg and in the Schrodinger picture rests ultimately on 
the following operator identity: 



-Y 



AY...AX 



where [F, . . .]X = Thus, in quantum-classical 

theory, one would like to derive an operator identity anal- 
ogous to that in Eq. (|13p . However, as already shown in 
Ref. [3| , because of the non associativity of the quantum- 
classical bracket in Eq. the identity that can be de- 
rived is 



= s 



e K 



(14) 



where the two operators 



n 



H 



= H- 



(15) 
(16) 



have been introduced and S is an ordering operator which 
is chosen so that the left and the right hand side of 
Eq. (HD) coincide by construction [3] , when the exponen- 
tial operators are substituted with their series expansion. 
The existence of such an ordering problem, and of the 
ordering operator S, in Eq. (fH|) is caused by the Pois- 
son bracket parts of the operators in Eqs. and (|16p . 
Hence, one can imagine that the solution to this problem 
can be found by dealing properly with these parts of the 
brackets. To this end, one can consider the quantum- 
classical equation of motion for the density matrix 



dp 

dt 



I r - 
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(17) 



As above discussed, in Eq. p?|) the ordering problem 
arises from the terms in the right hand side containing the 
Poisson bracket operator {...,.. .j^. Then, considering 
the identity 1 = P'P^^ =^ P^^'Pi Eq. (jl7p can be rewritten 
as 



dtp = - 



[H i] 


2i 



1 



'1 + I {/5, • ■ 



-^[H pp-n 




1 _ A I 



1 + |{P^---}. 
^ H 
P . 





' h' 




i 



H 

p-^p 



(18) 



where 
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^A---Mp)}, 





(19) 



(13) 

The operator X>e,[p] in Eq- CM depends from the 



quantum-classical density matrix, p, itself. However, if 
one momentarily disregards this non-linear dependence, 
Eq. ([18]) can be manipulated algebraically in order to de- 
velop a wave picture of quantum-classical mechanics. To 
this end, one can introduce quantum-classical wave fields, 
\ip{X)) and {ip{X)\, and make the following ansatzhr the 
density matrix 



p{X) = ^^J^XX))(^XX)| 



(20) 
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where one has assumed that, because of thermal disor- 
der, there can be many microscopic states \ip'-{X)) [i — 
which correspond to the same value of the 
macroscopic relevant observables In terms of the 

quantum-classical wave fields, \il}''{X)) and and 
considering the single state labeled by t, Eq. (llSp becomes 



W{x)){ij^{x)\ + |^^(x))(v>(x)| = 

+ W{X)){i^^{X)\H) 

+ i({ff,ln(p)}jV'(X))(^'(X)| 

- W{X)){t{X)\[Hp),H] 



(21) 

Equation (PT|) can be written as a system of two coupled 
equations for the wave fields [isj : 



^^J^\'<P(x,t)) 



H 



2i 



d 



X,t)\ 



(22) 



Equations (P^ . which are obeyed by the wave fields, are 
non-linear since their solution depends self-consistently 
from the density matrix defined in Eq. ([201) . These 
equations are also non-Hermitian since the operators 

|i?,ln(p)| and |ln(p), 7f I are not Hermitian. How- 
ever, this does not cause problems for the conservation 
of probability. The wave fields and evolve ac- 
cording to the different propagators 



^B,[p](0 = exp 
^sjp](i) = exp 



so that time-propagating wave fields are defined by 

\^P'{x,t)) = Z7,,[^](<)|v/(X)) 
(^^(x,i)| = (^^(x,)|l7,jp](t) . 

Quantum classical averages can be written as 

(x)(t) = [ dxJ2^^{^'ix,t)\x\rix,t)) 



(23) 
(24) 

(25) 
(26) 

(27) 



One can always transform back to the operator picture 
to show that the probability is conserved. 



Non-linear wave dynamics by means of 
non-Hamiltonian brackets 



The wave equations in ()22|) were derived starting from 
the non-Hamiltonian commutator expressing the dynam- 
ics of phase space dependent operators ■ It is interest- 
ing to recast quantum-classical wave dynamics itself by 
means of non-Hamiltonian brackets. It turns out that 
this form of the wave equations generalizes the mathe- 
matical formalism first proposed by Weinberg in or- 
der to study possible non-linear effects in quantum me- 
chanics (see Appendix |^ . 

Consider a case in which a single state is present, i.e. 
L = 1. Then, consider the wave fields \'4>) and (V'l as 
coordinates of an abstract space, and denote the point of 
such a space as 



c = 



Introduce the function 

n 



1^) 



(28) 



(29) 



and the antisymmetric matrix operator 



-1 





n {ln(p),g}„l^) 

2j (VliJ" 







Equations (|22|) can be written in compact form as 



(30) 



dc 

dt 



on 

dW) 



on 



dW') 



(31) 



Equations ((22|) . or their compact "Weinberg-like" form 
in Eq. pip , express the wave picture for the quantum- 
classical dynamics of phase space dependent quantum de- 
grees of freedom 0, Q . Such a wave picture makes one 
recognize the intrinsic non-linearity of quantum-classical 
dynamics. This specific features will be discussed, among 
other issues, in the next section. 



IV. ADIABATIC BASIS REPRESENTATION 
AND SURFACE-HOPPING SCHEMES 



Equations ([22)1 are written in an abstract form. In 
order to devise a numerical algorithm to solve them, one 
has to obtain a representation in some basis. Of course, 
any basis can be used but, since one would like to find a 
comparison with surface-hopping schemes, the adiabatic 
basis is a good choice. To this end, consider the following 
form of the quantum-classical Hamiltonian operator: 



H 



2M 



h{R) , 



(32) 
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where the first term provides the kinetic energy of the 
classical degrees of freedom with mass M, while h{R) 
describes the quantum sub-system and its coupling with 
the classical coordinates R. The adiabatic basis is then 
defined by the following eigenvalue equation: 

h\a; R) ^ Ea{R)\a; R) . (33) 
Since the non-linear wave equations in (j22p have been de- 



rived from the bracket equation for the quantum-classical 
density matrix (|17p , by dealing in a suitable manner with 
the Poisson bracket terms, the most simple way to find 
the representation of the wave equations in the adi- 
abatic basis is to first represent Eq. (fTT]) in such a basis 
and then deal with the terms arising from the Poisson 
brackets. The adiabatic representation of Eq. (fTT)) is 3] 

dtPaa'{X,t) = -^i£aa',l3p'Pl3l3'{X,t) , (34) 
00' 



where 



i^aa'.PP' — i^^aa\l3[j'^al3^a'l3' — Jaa',00' 

— {i^aa' + iLaa') SafjSa' f3' — Jaa',l3p' (35) 

Here, uoaa' = {Ea{R) — Ea'{R)) /h = Eaa'/^ and 



where 



F^^^{a:R\^^\a;R) 



OR 



(36) 



(37) 



is the Hellmann-Feynman force for state a. The operator 
J that describes nonadiabatic effects is 

P ^ / 1 ^ 9 \ , 

Jaa',l3fj' - - dap I 1 + -j^^af} ' ^ I Oa'p' 

(38) 

where dap — (a; R\{d / dR)\(3] R) is the nonadiabatic cou- 
pling vector and 



P 

Sap = Eapdap ( — ' dafi 



(39) 



Using Eqs. ([36|) and (f38|) . the equation of motion for 
the density matrix in the adiabatic basis can be written 
explicitly as 

M oR 



1 d 



E- 



d_ 



(40) 



The wave fields \'4)''{X)) and {'t]j''{X)\ can be expanded in 
the adiabatic basis as 

W'{X)) = ^|a;i?)(a;i?|V''(X))=^da;i?) 

Q a 

(V'^(X)I = ^(V'>;i?)(a;i?|=^(a;i?|Cr(X), 

a Q 

(41) 

and the density matrix in Eq. (I20p becomes 

(X, = E *)^^' ^) • (42) 

In order to find two separate equations for and C^*, 
one cannot insert Eq. ([42]) directly into Eq. (|40|) because 
of the presence of the derivatives with respect to the 
phase space coordinates R ad P. One must set Eq. ([40]) 
into the form of a multiplicative operator acting on paa' ■ 
To this end, for example, consider 

^P/3a' = E i^PP^ "^T"' " E i'^PM ^wVp"' 



E 



9(ln/5)/ jM 
dP 



7P 



'Plia' 



(43) 



Equation ([43]) shows how to transform formally a deriva- 
tive operator acting on p into a multiplicative operator 
which, however, depends on p itself. Therefore, Eq. ((40|) 
becomes 



"IT EaPaa' ^" ~rEa'Pa.a' 

n h 

^ p ^ p 



M 

1 _P d{\n p)a^ 
2^m' or 



1 p 5(lnp)a^ 



2 ^ Af aP 



'Pafj. 



2^^" dP 



a(lnp)^Q/ 



dP 



2^ M 



(44) 



Inserting the adiabatic expression for the density matrix, 
given in Eq. ((42|) . into Eq. (|44|l . one obtains, for each 
quantum state i, the following two coupled equations 

ciix,t) = -lij„c;,(x,t)-^^.d„^c;;(x,i) 
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2^ M 



1 ^ P a(inp)^ 



(45) 



2 ^ M 



2^M- dR ^^^^'^> 
^E^--^^^^r(^^0- (46) 



Quantum-classical averages of arbitrary observables can 
be calculated in the adiabatic as 



(x)W 



E 

0.0.' 



dXC'^{X,t)C'^*,{X,t)xc.'c.{X) , 



(47) 

where the coefficients {X, t) and C^* (X, t) are evolved 
according to Eqs. (|45p and (I46p . respectively. Equa- 
tions (|45p and (j46p are non-linear equations which couple 
all the adiabatic states used to analyze the system. 

At this stage, a general discussion about such a non- 
linear character is required. With a wide consensus, 
quantum mechanics is considered a linear theory. This 
leads, for example, to the visualization of quantum tran- 
sitions as instantaneous quantum jumps. The linearity 
of the theory also determines the need of considering in- 
finite perturbative series which must be re-summed in 
some way in order to extract meaningful predictions. 
Density Functional Theory is an example of a non-linear 
theory [19] but it is usually considered just as a com- 
putational tool. However, there are other approaches to 
quantum theory that represent interactions by an intrin- 
sic non-linear scheme f20]. It is not difficult to see how 
this is possible. Matter is represented by waves, these 
very same waves enter into the definition of the fields 
defining their interactioii [2l|. This point of view has 
been pursued by Jaynes and Barut [11], among oth- 
ers. These non-linear approaches depict quantum transi- 
tions as abrupt but continuous events [13] in which, to go 
from state ]1) to state ]2), the system is first brought by 
the interaction in a superposition a]l)-|-/3]2), and then, as 
the interaction ends, it finally goes to state ]2). It is un- 
derstood that this is made possible by the non-linearity of 
such theories because, instead, a linear theory would pre- 
serve the superposition indefinitely. Incidentally, the pic- 
ture of the transition process just depicted also emerges 
from the numerical implementation Q of the nonadia- 
batic quantum-classical dynamics of phase space depen- 
dent operators [3, [I]- The action of the operator J in 



Eq. ([35)) can build and destroy coherence in the system 
by creating and destroying superposition of states. As 
explained above, this is a feature of a non-linear theory. 
Such a non-linear character is simply hidden in the op- 
erator version of quantum-classical dynamics and clearly 
manifested by the wave picture of the quantum-classical 
evolution, which has been introduced in this paper. 

Since Eqs. ([IS]) and (|46p arc non-linear, their numeri- 
cal integration requires either to adopt an iterative self- 
consistent procedure (according to which one makes a 
ffi-st guess of paa', as dictated by Eq. (|42|l . calculates 
the evolved C''^{X, t) and C^*(X, t), and then goes into a 
recursive procedure until numerical convergence is ob- 
tained) or to choose a definite form for pj^^/, follow- 
ing physical intuition, and then calculating the time 
evolution, according to the form of Eqs. (|l5)) and ([^5]) 
which is obtained by using p^^/. This last method is 
already known within the Wigner formulation of quan- 
tum mechanics ,24] as the method of Wigner trajecto- 
ries (25| . It is also important to find some importance 
sampling scheme for the phase space integral in Eq. (j47p . 
Such sampling scheme may depend on the specific form 
Xaa' of the observable. It is interesting to note that 
Eqs. (gSl), (gll), and (gT]) can be used to address both 
equilibrium and non-equilibrium problems on the same 
footing. However, the dynamical picture provided by 
Eqs. (US]) and is very different both from that of 
the usual surface-hopping schemes [2^ and from that of 
the nonadiabatic evolution of quantum-classical opera- 
tors In order to appreciate this, for simplicity, one 
can consider a situation in which there is no thermal dis- 
order in the quantum degrees of freedom so that u — 1: 
viz., the density matrix becomes that of a pure state 
Paa'{X,t) Cc{X,t)Cl,{X,t). Then, equations gS]) 
and (I46p remain unaltered and one has just to remove the 
index l from the coefficients. Therefore, it can be real- 
ized that no classical trajectory propagation, and no state 
switching are involved by Eqs. ([^5]) and (|46p . Instead, in 
order to calculate averages according to Eq. ([Tf)) . one has 
to sample phase space points and integrate the matrix 
equations. 

In the next section, an equilibrium approximation of 
Eqs. ([45l) and ([46]), along the lines followed by the method 
of Wigner trajectories |25i], is given and applied, with 
good numerical results, to the adiabatic and nonadiabatic 
dynamics of the spin-boson model. 



V. WAVE DYNAMICS OF THE SPIN-BOSON 
MODEL 

The theory developed in the previous sections can be 
applied to simulate the relaxation dynamics of the spin- 
boson system [27|. This system has already been stud- 
ied within the framework of quantum-classical dynamics 
of operators in Ref. [13] and "exact" numerical results 
were obtained at short-time by means of an iterative 
path integral procedure developed by Nancy Makri and 
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co-worker ^2§\ . The short-time results of Ref. [27*1 numer- 
ically coincide with those obtained by the path integral 
calculation of Ref. However, as it is shown later 

by Fig. [21 the quantum-classical results of Ref. ^2?| have 
some limitations concerning the numerical stability of the 
algorithm beyong a certain time length. Using the di- 
mensionless variables of Ref. fl^, the quantum-classical 
Hamiltonian operator of the spin-boson system reads 



The process of relaxation from the initial state can be fol- 
lowed by monitoring the subsystem observables dz , which 
in the adiabatic basis reads 



1 ( 2G 1 - G2 



(56) 



The adiabatic basis is real so that the initial density 
matrix of the system can be written as 

2 

p„„.(X,0) = ^7A„(X,0),^„.(X,0), (57) 



hs + Hb + Vc{R) , 



(48) 



where hs = — iliTi is the subsystem Hamiltonian, Hi, — 
EU ^'/2 + l/2c.2i?2 = J2f=i + V,{R) is the 

Hamiltonian of a classical bath of N harmonic oscilla- 
tors, and Vc{R) = — '^j^zRj = l{R)dz is the inter- 
action between the subsystem and the bath. An Ohmic 
spectral density is assumed for the bath. Hence, denot- 
ing the Kondo parameter as S^k and the cut-off frequency 
as Wmaxi the frequencies of the oscillators are defined by 
ujj = -ln(l - jwo), where = iV"^(l - exp(-Wmax)), 
and the constants entering the coupling by and Cj = 
\/£,K^^o ^j- The adiabatic eigenvalues and eigenvectors, 
respectively, are 



Ei,2 ^VbT ^n^+i^{R) 



(49) 



|2;i?) = 



1 / 1 + G 
v/2(l + G2 \ l-G 

1 / -1 + G 



where 



G{R)^^-\R) -n + + -f^{R) 



(50) 



(51) 



The coupling vector daa' — (a; ^| d /dR\a'; R) is 
di2=-d2i = (l + G^y^dG/dR. 



Assuming an initially uncorrelated density matrix, where 
the bath is in thermal equilibrium and the subsystem is 
in state | t), the initial quantum-classical density matrix 
in the adiabatic basis takes the form 



where 
PsiO) 

and 



p(o) = pMPbix) , 

1 f{l + G)^ 1-G2 
2(1 + G2) [ 1~G^ (1-G)2 



(53) 
(54) 



N 



pbix) = n 



1=1 

X exp 



tanh(/3wj/2) 

2tanh(/3w,/2) f if Jjl^^^ 



(55) 



where 



V'i(x,o) = ^i(x,o) = 



1 + G 



V'2(X,0) =02(X,O) - ^Pb{X) 



v/2(l + G2) ' 

1-G 
V2(l + G2) 



(58) 



(59) 



Such coefficients enter into the calculation of the observ- 
able 

{(Tz{t)) = y rf^0a'(X,tX'"(X)^,(X,t) . (60) 



The coefficients evolve in time according to Eqs. ps)) 
and (|46)) . where one must set G^ = tpa and C-*, = (fia' ■ 
In order to devise an effective computational scheme for 
such equations, one could assume that the density ma- 
trix entering Eqs. (^5)) and is taken to be that at 
t = oo, when the total system (subsystem plus bath) has 
reached thermal equilibrium. The equilibrium quantum- 
classical density matrix is known as a series expansion in 
If one makes the additional assumption of complete 
decoherence at t = 00, only the 0{h'^) term can be taken 



,(0) 



(52) where Zo = Eaa' / dX pf (X). Then 



9 In p, 



{0)aa' 



dR 



-fi-Q^5^c.'=PF^{R)5c.c'. (62) 



9 In p, 



(O)Qa' 



dP 



(63) 



Equations ([l5)) and become 
d 



dt 



Mx,t) = -iEMX,t) 



J2P-dcpil-^E^0)MX,t) (64) 



-6„,{X,t) = iE^,(j)^,{X,t) 
dt 



(65) 
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In Eqs. (gel) and §Q the terms ±(/3/2)P • (and 
the analogous terms with S,a') cancel each other. In the 
adiabatic basis = d22{R) — 0. Hence, defining 

the matrix 



-iEi 

P ■ di2 (l + I-B12 



-P • dl2 ( 1 

-iE2 



(1 - f £^12) 



(66) 



Equations (|64p and ((65)) can be written as 



d_ 

di 



■01 
1p2 



= s • 



"02 



(67) 



which can be integrated by means of the simple algo- 
rithm *(X, dr) = *(X, 0) + dT&{X, 0) • qiX, 0), where 
'3/ = ("0, 0) and = (S, S*). The phase space part of 
the initial values of and (f) can be used as the weight 
for sampling the coordinates X entering the classical inte- 
gral in Eq. (pD|) . Then, for each initial value X, Eqs. ([^ 
must be integrated in time so that averages can be cal- 
culated. It is worth to note that in such a wave scheme 
the Eulerian point of view of quantum-classical dynam- 
ics [1, is preserved. This is different from what hap- 
pens in the original operator approach @, , where in 
order to devise an effective time integration scheme by 
means of the Dyson expansion, one is forced to change 
from the Eulerian point of view (according to which the 
phase space point is fixed and the quantum degrees of 
freedom evolve in time at this fixed phase space point) 
to the Lagrangian point of view, where phase space tra- 
jectories are generated. Moreover, it must be noted that 
the numerical integration of Eqs. (|67p provides directly 
the nonadiabatic dynamics without the need to introduce 
surface-hopping approximations. 

In order to be able of comparing the results with those 
presented in Ref . [To| , the numerical values of the param- 
eters specifying the spin-boson system have been chosen 
to be /3 = 0.3, n ^ 1/3, w,nax ^ 3, ^ 0.007, and 
A'' = 200. Figure [T] shows the results in the adiabatic 
case, obtained by setting du = in Eqs. (|57)) . One can 
see that, in spite of the simple approximation of the form 
of the density matrix made in the equations of motion, 
the wave theory provides results which are in good agree- 
ment with those obtained with the operator approach of 
Ref. Instead, Fig.[2]shows the results of the nonadia- 
batic calculation. This is to be compared with the results 
of the operator theory [lol (which are identical with the 
exact" ones of Ref. [23 )■ Of course, since different ways 
of dealing with the nonadiabatic effects are used in the 
two approaches the results do not need to be the same. 
However, the results of the wave theory follow qualita- 
tively those of Ref. while improving substantially the 
statistical convergence and increasing the length of the 
time interval spanned by a factor of 2— 3. Such results are 
particularly encorauging and suggest the possible appli- 
cation of the wave theory here proposed, for example, to 



0.5 



-0.5 




FIG. 1: Adiabatic dynamics of the spin-boson model. f3 = 0.3, 
n = 1/3, LJmax = 3, Cif = 0.007, iV = 200. The black circles 
show the results of the calculation with the theory proposed 
in this paper while the light dashed line with squares shows, 
for comparison, the results of the calculation in Ref. [lol |. 




FIG. 2: Nonadiabatic dynamics of the spin-boson model. 
f3 = 0.3, n = 1/3, tJmax = 3, = 0.007, iV = 200. The black 
circles show the results of the calculation with the theory pro- 
posed in this paper while the line with error bars shows, for 
comparison, the results of the calculation in Ref. [l3|- These 
latter are indistinguishable, at short time, from the "exact" 
results of Ref. [28l ]. which were obtained by means of an iter- 
ative path integral procedure. 



the calculation of nonadiabatic rate constants of complex 
systems in the condensed phase 



VI. CONCLUSIONS 

In this paper the approach to the quantum-classical 
mechanics of phase space dependent operators has been 
remodeled as a non-linear formalism for wave fields. It 
has been shown that two coupled non-linear equations for 
phase space dependent wave fields correspond to the sin- 
gle equation for the quantum-classical density matrix in 
the operator scheme of motion. The equations of motion 
for the wave fields have been re-expressed by means of a 
suitable bracket and it has been shown that the emerging 
formalism generalizes within a non-Hamiltonian frame- 
work the non-linear quantum mechanical formalism that 
has been proposed recently by Weinberg. Finally, the 
non-linear wave equations have been represented into the 
adiabatic basis and have been applied, after a suitable 
equilibrium approximation, to the numerical study of the 
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adiabatic and nonadiabatic dynamics of the spin-boson 
model. Good results have been obtained. In particular, 
the time interval that can be spanned by the nonadiabatic 
calculation within the wave scheme of motion turns out 
to be a factor of two-three longer than that accessible 
within the operator scheme of motion. This encourages 
one to pursue the application of the wave scheme of mo- 
tion to the calculation of correlation functions for systems 
in the condensed phase. Future works will be specifically 
devoted to such an issue. 
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dV ' 



i dH 

hd\¥) 



Bi 



(A6) 



It is easy to see that, when the Hamiltonian function is 
chosen as in Eq. jASj, Eq. (jA4|, or its explicit form (jA5I 
IA6[) . gives the usual formalism of quantum mechanics. 
It is worth to remark that in order not to alter gauge 
invariance, the Hamiltonian and the other observables 
must obey the homogeneity condition: 



n = {mdn/d(:2)) = {{dH/dCim 



(A7) 



Weinberg showed how the formalism above sketched can 
be generalized in order to describe non-linear effects in 
quantum mechanics [T^ . To this end, one must maintain 
the homogeneity condition, Eq. (jA7p , on the Hamiltonian 
but relax the constraint which assumes that the Hamilto- 
nian must be a bilinear function of the wave fields. Thus, 
the Hamiltonian can be a general function given by 



(A8) 



APPENDIX A: WEINBERG'S FORMALISM 

Consider a quantum system in a state described by 
the wave fields j^*) and (^'|, where Dirac's bra-ket nota- 
tion is used to denote *(r) = {r\^i) and **(r) = (*|r). 
Observables are defined by functions of the type 



a=(^\A\^) , 



(Al) 



where the operators are Hermitian, A = . Weinberg's 
formalism can be introduced by defining Poisson brackets 
in terms of the wave fields j^*) and To this end, one 
considers the wave fields as "phase space" coordinates 
C = d^"), (*|), so that Ci = |*> and C2 = (^'1, and then 
introduce brackets of observables as 



{a, 6}e 



Eda_ db_ 



(A2) 



The bracket in Eq. \K2\ defines a Lie algebra and a 
Hamiltonian systems. Typically, the Jacobi relation 
is satisfied, i.e. J = {a, {&, cjelg + {c{a,b}^}^ + 
{6, {c, a}g}g = 0. In order to obtain the usual quantum 
formalism, one can introduce the Hamiltonian functional 
in the form 



(A3) 



where H is the Hamiltonian operator of the system. 
Equations of motion for the wave fields can be written in 
compact form as 



(A4) 



The compact form of Eq. (jA4p can be set into an explicit 

form as 



i &H 



rid{^\ 



B2 



(A5) 



where n is arbitrary integer that fixes the order of the 
correction, Tio = ^, and 

Hi = j drdr'dr"dr"'-^*{r)^*{r') 

X G(r,r',r",r"')*(r")*(r"') , (A9) 

with analogous expressions for higher order terms. Appli- 
cations and thorough discussions of the above formalism 
can be found in Ref. [T^ . 

Once Weinberg's formalism is expressed by means of 
the symplectic form in Eq. (|A4[) . it can be generalized 
very easily in order to obtain a non-Hamiltonian quan- 
tum algebra. To this end, one can substitute the an- 
tisymmetric matrix B with another antisymmetric ma- 
trix CI = ^[C]i whose elements might be functionals 
of C = (l^*): (^'D obeying the homogeneity condition in 
Eq. (|A7p . By means of Jl a non-Hamiltonian bracket 
{...,.. .}^ can be defined as 



r ,1 (^o. ^ db 



(AlO) 



In general, the bracket in Eq. (|A10[) does no longer satisfy 
the Jacobi relation 

J = {a, {b, c}J^ + {c {a, b}J^ + {b, {c, a}J^ ^ . 

(All) 

Thus, non-Hamiltonian equations of motion can be writ- 
ten as 



dt ~ h 



(A12) 



In principle, the non-Hamiltonian theory, specified by 
Eqs. (jAlOp . (jAlip . and (|A12p . can be used to address the 
problem of non-linear correction to quantum mechanics. 
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as it was done in Refs. [l^- In the present paper, it has 
been shown that such a non-Hamiltonian and non-hnear 
version of quantum mechanics is already imphed when 
one formulates quantum-classical dynamics of operators 



by means of suitable brackets. As a matter of fact, it was 
shown that the quantum-classical theories of Refs. 0, Q 
can be mapped onto a wave formalism which has precisely 
the same form specified by Eqs. (|XTO)) . (|XTT|) . and ([Xl2|) . 
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